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Abstract 

The measurement of sin2/3 by the BABAR experiment, where (3 is one of the 
angles of the Unitarity Triangle, is described. Some prospects for the future 

of the measurement are also discussed. 
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1 Introduction 



The BABAR experiment consists of an asymmetric electron-positron collider operating at 
the T(4S) resonance. More details on the detector can be found in [||. The aim is to 
overconstrain the unitarity triangle by measuring its sides and angles. The analysis reviewed 
here measures sin2/3 by studying time-dependent CP violating asymmetries in B° — > J/ipK® 
and B° — > ip'K° s decays. 

2 Overview of the sin2/3 analysis 

There are five main parts to measuring the CP violating asymmetry: 

• Selection of signal CP events 

• Measurement of the distance Az between the two B° decay vertices along the T(4S) 
boost axis 

• Determination of the flavour of the tag-side B 

• Measurement of dilution factors for the different tagging categories 

• Extraction of sin2/? via an unbinned maximum likelihood fit 

2.1 Event Selection 

The sample used for the analysis is 9.8 fb" 1 of data recorded between January and July 2000 
of which 0.8 fb _1 was recorded 40 MeV below the T(4S) resonance. Particle identification 
uses mainly the Csl calorimeter for electrons, the Instrumented Flux Return for muons and 
the DIRC for kaons. Extra information is provided by dE/dx measured in the tracking 
system. The selection for the CP events proceeds as follows. Pairs of electrons or muons 
coming from a common vertex are combined to form J/ip and ip' candidates. The ip' is a ls° 
reconstructed from its decay into J/ip -k + -k^ . The K$ candidates are made from either a 
pair of charged tracks or a pair of 7r° candidates. In addition there are various event shape 
and topological cuts designed to reduce continuum and BB background. Full details of the 
selection can be found in 0. The final event sample is shown in figure [3]. 

There are two other B decay samples. One consists of fully reconstructed semileptonic 
(B° — > D*~l + vi) and hadronic (B° — > D^~ir + , D^~p + , D^~af) decays as well as a control 

(*)0 i 

sample of B + — > D ir + events. The selection of this sample is described in g] and . The 
other is a charmonium control sample containing fully reconstructed neutral or charged B 
candidates in two-body decay modes with a J/ip in the final state (e.g. B + — > J/ip K + , B° — > 
J/ip K*°(K*° -> K + tt-)). 
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Figure 1: CP signal event distributions for J/il)Ks(7t + 7r ) (left), J/ipKs(iT TT ) (middle) and 
iff' K s (7T + ir~) (right). 

2.2 Measuring Az 

The time-dependent decay rate for the E>cp is given by 

/±(r, Am d , V sin 2f3, t) = ^re" r|f| [1 ± I? sin 2/? x sinAm d t] (1) 

where the + or - sign indicates whether the B tag was tagged as a B° or B° respectively The 
dilution factor T> is given by T> = 1 — 2w, where w is the mistag fraction (the probability 
that the B tag is identified incorrectly). To account for finite detector resolution, the time 
distribution must be convoluted with a resolution function: 



K(Az;a) = £ _^ = exp(-(A^-^) 2 /2^ 2 ) , (2) 

i=1 CW27T V ' 

which is just the sum of two Gaussians where the fc, Si and Oi are the normalizations, biases 
and widths of the distributions. In practice two scale factors S\ and S2 are introduced such 
that <7j = Si x a^t where a^t is an event-by-event calculated error on At. They take account 
of underestimating the uncertainty on At due to effects such as hard scattering and possible 
underestimation of the amount of material traversed by the particles. The resolution function 
parameters are obtained from a maximum likelihood fit to the hadronic B° sample and are 
shown in table [I]. The f w parameter represents the width of a third Gaussian component, 
included to accommodate a small (~1%) fraction of events which have very large values of 
Az, mostly caused by vertex reconstruction problems. This Gaussian is unbiased with a 
fixed width of 8 ps. Further details can be found in . 

2.3 B flavour tagging 

Each event with a CP candidate is assigned a B° or B° tag if it satisfies the criteria for 
one of the several tagging categories. The figure of merit for each tagging category is the 



Parameter 


Value 


<*i (ps) 


-0.20 ±0.06 


from fit 


Si 


1.33 ±0.14 


from fit 


U (%) 


1.6 ±0.6 


from fit 


h (%) 


75 


fixed 


S 2 (ps) 





fixed 


s 2 


2.1 


fixed 



Table 1: Resolution function parameters. Those, labeled 'from fit' are measured from data 
and those marked 'fixed' are determined from Monte Carlo. 

effective tagging efficiency Qi = ej(l — 2wi) 2 where is the fraction of events assigned to 
category i and Wi is the probability of mis-tagging an event in this category. The statistical 
error on sin2/9 is proportional to 1/yfQ where Q = Y^iQi- There are five tagging categories: 
Electron, Muon, Kaon, NT1 and NT2. 

The first three require the presence of a fast lepton and/or one or more charged kaons 
in the event and depend on the correlation between the charge of a primary lepton or kaon 
and the flavour of the b quark. If an event is not assigned to either the Electron or Muon 
categories, it is assigned to the Kaon category if the sum of the charges of all the identified 
kaons in the event is different from zero. If both lepton and kaon tags are available but 
inconsistent the event is rejected from both categories. 

NT1 and NT2 are categories from a neural network algorithm, this approach being mo- 
tivated by the potential flavour-tagging power carried by non-identified leptons and kaons, 
correlations between leptons and kaons and more generally the momentum spectrum of 
charged particles from B meson decays. The output of the neural network tagger xnt can 
be mapped onto the interval [-1,1] with xnt < representing a B° tag and xnt > a B° tag. 
Events with \xnt\ > 0.5 are classified in the NT1 category and events with 0.2 < \x^t\ < 0.5 
in the NT2 category. Events with \x NT \ < 0.2 are excluded from the final analysis sample. 

2.4 Measurement of tagging performance 

The effective tagging efficiencies and mistag fractions for all the categories are measured 
from data using a maximum likelihood fit to the time distributions of the B° hadronic 
event sample. The procedure uses events which have one B fully reconstructed in a flavour 
eigenstate mode. The tagging algorithms are then applied to the rest of the event, which 
represents the potential B tag - Events are classified as mixed or unmixed depending on whether 
the B tag is tagged with the same or opposite flavour as the Bcp. One can express the 
time-integrated fraction of mixed events x as a function of the B°B° mixing probability, 
X — Xd + (1 — ^Xd)w where Xd = 1^1/(1 + x d)i w ^ n x d — Am d /r. Thus an experimental 
value of the mistag fraction w can be deduced from the data. 

A more accurate estimate of w comes from a time-dependent analysis of the fraction of 
mixed events. The mixing probability is smallest at low At so that this region is governed 



by the mistag fraction. Figure |2| shows the fraction of mixed events versus At. The resultant 
tagging performances are shown in table |2|. 
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Figure 2: The fraction of mixed events as a function of |At| for data events in the hadronic 
sample for neutral B mesons (full squares) and charged B mesons (open circles). The dot- 
dashed line at t cut = 2.5 ps indicates the bin boundary for the time-integrated single-bin 
method. 



2.5 Extracting sin2/3 

A blind analysis technique was adopted for the extraction of sin2/3 to eliminate possible 
experimenter bias. The technique hides both the result of the likelihood fit and the visual CP 
asymmetry in the At distribution. This method allows systematic studies to be performed 
while keeping the numerical value of sin2/3 hidden. 

Possible systematic effects due to uncertainty in the input parameters to the fit, incom- 
plete knowledge of the time resolution function, uncertainties in the mistag fractions and 
possible limitations in the analysis procedure were all studied. Details can be found in [[|. 
The systematic errors are summarized in table [5| 

2.6 Results and checks 

The maximum likelihood fit for sin2/?, using the full tagged sample of 120 B° — > J/tpK® and 
B° — > ip'K® events yields: 



Tagging category 


e(%) 


w (%) 


Q(%) 


Lepton 


11.2 ±0.5 


9.6 ± 1.7 ± 1.3 


7.3 ±0.3 


Kaon 


36.7 ±0.9 


19.7 ± 1.3 ± 1.1 


13.5 ±0.3 


NT1 


11.7 ±0.5 


16.7 ±2.2 ±2.0 


5.2 ±0.2 


NT2 


16.6 ±0.6 


33.1 ±2.1 ±2.1 


1.9 ±0.1 


all 


76.7 ±0.5 




27.9 ±0.5 



Table 2: Tagging performance as measured from data. 



Source of uncertainty 


Uncertainty on sin2/9 


uncertainty on 


0.002 


uncertainty on Am d 


0.015 


uncertainty on Az resolution for CP sample 


0.019 


uncertainty on time-resolution bias for CP sample 


0.047 


uncertainty on measurement of mistag fractions 


0.053 


different mistag fractions for CP and non-CP samples 


0.050 


different mistag fractions for B° and B° 


0.005 


background in CP sample 


0.015 


total systematic error 


0.091 



Table 3: Summary of systematic uncertainties. The different contributions are added in 
quadrature. 
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Figure 3: Variation of the log likelihood as a function of sin2/3. The two horizontal dashed 
lines indicate changes in the log-likelihood corresponding to one and two statistical standard 
deviations. 



sin2/3 = 0. 12 ± 0.37 (stat) ± 0.09 (syst) (preliminary) . (3) 

The log likelihood is shown as a function of sin2/3 in figure |3|. The raw asymmetry as a 
function of At is shown in figure [| 

The probability of obtaining a statistical uncertainty of 0.37 is estimated by generating a 
large number of toy Monte Carlo experiments with the same number of tagged CP events as 
in the data sample. The errors are distributed around 0.32 with a standard deviation of 0.03, 
meaning that the probability of obtaining a larger statistical error that the one observed is 
5%. From a large number of full Monte Carlo simulated experiments, we estimate that the 
probability of finding a lower value of the likelihood than the one observed is 20%. 

Several cross-checks are performed to validate the main analysis. The charmonium and 
fully-reconstructed hadronic control samples are composed of events that should exhibit no 
time-dependent asymmetry. These events are fitted in the same way as the signal CP events 
to extract an "apparent CP asymmetry" . The results are shown in table 01. 



2.7 Constraints on the Unitarity Triangle 



The Unitarity Triangle in the (p, 77) plane is shown in figure [5[ The two solutions cor- 
responding to the measured central value are shown as straight lines. The cross-hatched 
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Figure 4: The raw B° - B° asymmetry (N B o - N So )/(N B o + N Eo ). The time-dependent 
asymmetry is represented by a solid curve for the central value of sin2/9, and by two dotted 
curves for the values at plus and minus one statistical standard deviation from the central 
value. The curves are not centered at (0,0) because the CP sample contains an unequal 
number of B° and B° events (70 B° versus 50 B°). The \ 2 between the binned asymmetry 
and the result of the maximum likelihood fit is 9.2 for 7 degrees of freedom. 



Sample 


Apparent CP asymmetry 


Hadronic charged B decays 


0.03 ±0.07 


Hadronic neutral B decays 


-0.01 ±0.08 


J/ipK + 


0.13 ±0.14 


J/iIjK*°(K*° -> K+7T-) 


0.49 ±0.26 



Table 4: Summary of systematic uncertainties. The different contributions are added in 
quadrature. 




Figure 5: Present constraints on the position of the apex of the Unitarity Triangle with the 
BABAR result indicated by the cross-hatched regions. 



regions correspond to one and two times the one-standard-deviation experimental uncer- 
tainty. The ellipses represent regions allowed by all other measurements that constrain the 
triangle. They are shown for a variety of choices of theoretical parameters. More details can 
be found in || . 

3 Future prospects 

The preceding pages describe only a preliminary measurement of sin2/5 by the BABAR exper- 
iment. More data will allow extra channels to be included in the final fit as well as providing 
more events for the currently used decay modes. The new channels will bring extra exper- 
imental and theoretical challenges with them. Such present and future issues are discussed 
in the sections that follow. 

3.1 Available Modes 

The B decay modes that have been used to measure sin2/3 up to now are clean in that 
they are vector-scalar, b —>■ ccs transitions which have no significant pollution from penguin 
diagrams. The next step is to add vector- vector modes such as B° — > J/tpK*. These modes 



require an angular analysis of the vector meson decay products, due to the different partial 
waves and therefore admixture of CP odd and CP even that is present in the final state. 
Such an angular analysis has already yielded preliminary results for the J/ipK* modes. Once 
one has measured the polarizations in these modes, they are as clean, theoretically, as the 
vector-scalar modes. Another obvious addition is B° — > J/ipK^ decays where the challenge 
here is to understand the background well enough to make the channel feasible. Work is 
ongoing in this area. 

A different kind of difficulty is presented by channels with a significant degree of penguin 
contamination, such as b — > ccd scalar-scalar modes (e.g. B° — > D + D~). Here the fit must 
take into account the fact that the true value of sin2/3 is shifted by an amount proportional 
to the ratio of tree to penguin contributions. This ratio is model dependent and subject to 
large theoretical uncertainties. 

Finally, modes such as B° — > D*D* and B — > J/ipp° which are vector- vector, b — > ccd 
transitions face the theoretical challenges of the penguin contaminated modes described 
above, as well as requiring an angular analysis to solve the vector-vector CP admixture 
problem. 

3.2 Experimental Considerations 

There are also experimental analysis issues which need to be resolved or studied in greater 
depth in the future. The tagging algorithms that BABAR uses should be developed and ex- 
tended to include extra tagging categories such as the using the soft pion from D* decays 
and incorporating leptons at an intermediate momentum (i.e. from a cascade). It would be 
useful to take account of correlations within an event, such as when two different tagging 
categories report an answer. This can give more information about the event if the correla- 
tions are well understood. There is also an open question when it comes to measuring the 
tagging performance from the hadronic or semileptonic B decay samples. One then needs 
to be absolutely sure that using exactly the same numbers for the CP signal event sample 
is a valid thing to do. 

The measurement of Az is another crucial part of the analysis and it is important that 
the errors and biases to this distribution are understood. The distribution tends to be biased 
by the decays of particles which fly significantly from the original B decay vertex, such as 
D°s. These can be rejected by looking explicitly for cascade decays. The parameterization 
of the resolution function incorporates detector effects such as misalignments and electronics 
readout effects. All contributions to the width should be studied in order to fully understand 
the error on Az. 

Backgrounds to the various CP modes can also be a problem. The channels vary in 
terms of how much background they experience and this background can be particularly 
dangerous if it has a significant structure in Az. For charmonium channels, much of the 
background comes from events containing a real J/ip. In that case, one needs to study 
exactly which modes contribute and what their shape is in the final distributions (if they 
cannot be removed otherwise). Non- resonant backgrounds to vector- vector modes such as 
the J/tl)K°ir° contribution to J/tl)K*°(K*° — > K°tt°) are in principle dangerous since they 



can have CP violating properties but no angular structure. However, the branching ratios 
for these non-resonant modes are typically poorly known and consistent with zero making 
it difficult to simulate them in the correct proportions. 

3.3 Study of Statistical Error 

It seems anomalous that both BABAR and Belle record higher statistical errors than one 
would expect. The fitting procedure is, and continues to be a vigorously studied part of the 
analysis as we need to be certain that the likelihood function is of exactly the correct form 
for the final fit. 

4 Conclusions 

A preliminary measurement of sin2/9 by BABAR has been presented. The errors on the final 
result make it difficult to express its significance in terms of constraints on the Unitarity 
Triangle. However, results based on a much larger data sample (~20 fb^ 1 ) will soon be 
available. Combined with a better understanding of systematic effects, this should make the 
next measurement of sin2/3 even more interesting than the current one. It is also expected 
that other CP modes will soon be available for analysis including B>° — > J '/t/jK*°(K*° — > 
K s tt°) and B° — > J/vpK L . The larger data sample with additional CP modes should yield 
a value of sin2/3 for which the statistical and systematic errors are about one-half of their 
current values. 
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